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Abstract 

We study damping of a dipole oscillation in a Bose-Condensed gas in a combined cigar-shaped 
harmonic trap and one-dimensional (ID) optical lattice potential at finite temperatures. In order 
to include the effect of thermal excitations in the radial direction, we derive a quasi-lD model of 
the Gross-Pitaeavskii equation and the Bogoliubov equations. We use the Popov approximation to 
calculate the temperature dependence of the condensate fraction with varying lattice depth. We 
then calculate the Landau damping rate of a dipole oscillation as a function of the lattice depth 
and temperature. The damping rate increases with increasing lattice depth, which is consistent 
with experimental observations. The magnitude of the damping rate is in reasonable agreement 
with experimental data. We also find that the damping rate has a strong temperature dependence, 
showing a sharp increase with increasing temperature. Finally, we emphasize the importance of 
the radial thermal excitations in both equilibrium properties and the Landau damping. 
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I. INTRODUCTION 



Recently the dynamics of ultracold atomic gases in optical lattices have attracted atten- 
tion both theoretically and experimentally [l]. In particular, center-of-mass dipole oscilla- 
tions of Bose-Einstein condensates in a combined cigar shaped trap and one dimensional 
(ID) optical lattice potential have been experimentally studied in detail In the 

presence of the periodic lattice potential, a decrease of the dipole mode frequency was ob- 
served {3)]. This decrease can be understood in terms of the increase of the effective mass 
due to the lattice potential 5j. On the other hand, at finite temperatures where an ap- 
preciable number of atoms are thermally excited out of the condensate, strong damping of 
;he dipole oscillation of the condensate was observed in the presence of the lattice potential 
4|. In a pure harmonic potential, a Bose gas exhibits undamped dipole oscillations even at 
finite temperatures, since the condensate and noncondensate atoms oscillate with the same 
frequency without changing their density profiles 6]. In contrast, in the presence of the 
periodic lattice potential, only the condensate component can coherently tunnel through 
the potential barriers, while the thermal component is locked by the lattice potential {4]. It 
is clear that this incoherent thermal conponent gives rise to the damping of the condensate 
oscillations. In Bose-condensed gases trapped in harmonic potentials, Landau damping is 
known to be the dominant contribution to damping of the condensate oscillations in the 
co.Msion.ess regime BQfl the analvs.s of damping of d.po.e oscillations in opti- 

cal lattice, the authors of Ref. [4( gave a rough estimate of the Landau damping rate and 
compared it with their experimental data. However, quantitative calculations of the Landau 
damping rate of dipole condensate oscillations in optical lattices have not given in any detail 
so far. In fact, even equilibrium properties of a Bose gas in an optical lattice at finite tem- 
peratures in connection with the experimentas of Ref. J4] have not been studied in detail. 
Although several papers have discussed finite-temperature properties of ultracold atoms in a 
ID optical lattice, most theoretical studies have concentrated on the first Bloch band using 
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This 



the Bose- Hubbard model, and have ignored the effect of radial excitations 11, 
approximation is only valid when k-^T -C -Er and k-^T <C huj_, where £r is the recoil energy 
that is roughly the highest energy in the first Bloch band, and hu± is the first excitation 
energy in the radial direction. In order to obtain more quantitative results that is applicable 
to the experiments of Refs. 3, J], however, it is important to consider thermal excitations in 
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the higher bloch bands as well asa the radial direction, since the experiments do not always 
satisfy k^T <C E-& and k-^T <C hu±. 

In this paper, we study the Landau damping of condensate oscillations of a trapped Bose 
gas in a ID optical lattice, with explicitly including the effect of the radial excitations. In 
Sec. II, we derive a quasi-lD model of the Gross-Pitaeavskii equation and the Bogoliubov 
equations that include the effect of the excitations in the radial direction. Using the Hatree- 
Fock-Bogoliubov-Popov (HFB-Popov) approximation, we solve these equations to calculate 
the temperature dependence of the condensate fraction. 

In Sec. Ill, we calculate the Landau damping of the dipole condensate oscillation in an 
optical lattice. We calculate the damping rate as a function of the lattice depth with a fixed 
temperature, and compare it with the experimental data 4]. The magnitude of the damping 
rate is found to be in reasonable agreement of the experimental data [4] . The damping rate 
increases with increasing lattice depth, which is consistent with the experimental result J4]. 
We also calculate the temperature dependence of damping rate with a fixed lattice depth. 
We find that the damping has a strong temperature, showing a sharp increase with increasing 
temperature. 



II. QUASI ID MODELING OF A TRAPPED BOSE GAS 



We consider a Bose condensed gas in a combined potential of highly-elongated harmonic 
trap and ID optical lattice. Our system is described by the following Hamiltonian: 

I 2 
2m 



H = J driver) 



-^V 2 + Kxt( 



+|^(r)^(r)^(r)^(r)|, (1) 

where the coupling constant determined by the s-wave scattering length 

a. The external potential Kxt is given by V ext (r) = V trap (r) + V op (z), where K rap (r) = 
y [o;5.(x 2 + y 2 ) + uj 2 z 2 } is an anisotropic harmonic potential and V op (z) = sE R cos 2 (kz) is an 
optical lattice potential. Here s is a dimensionless parameter describing the intensity of the 
laser beam creating the ID lattice in units of the recoil energy Er = where k = ^ is 
fixed by the wavelength A of the laser beam. 

In this paper we consider a highly-anisotropic cigar shaped harmonic trap potential u± ^> 
uj z . In order to take into account this quasi-lD situation, we expand the field operator in 
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terms of the radial wave function 14|: 

^(r) = E ipa(z)(j) a (x, y), (2) 

a 

where 4> a (x,y) is the eigenfunction of the radial part of the single-particle Hamiltonian, 

rr 



<Pa(x,y) 



(3) 



which satisfy the orthonormality condition J dxdy 0* (x, y) = 5 a p. Here a = (n x ,n y ) 

is the index of the single-particle state with the eigenvalue €( nx> n y ) — hu±(n x + n y + 1). 
Inserting Eq. (j2J) into Eq. ([Tj) and using Eq. we obtain 

ft 2 <9 2 



£ = E / 



rn 



+ E 



9aa'(3f3' 



aa'pP' 



(4) 



where the renormalized coupling constant is defined by 



WW = 9 / dxdy<f>* a (t>p<f>i3'(t>c 



(5) 



Following the procedure described in Ref. [15] , we separate out the condensate wavefunc- 
tion from the field operator as ip a = (ip^+ipa = &a + ipa, where $ Q = (ip a ) is the condensate 
wavefunction and ip a is the noncondensate field operator. Within the HFB-Popov approxi- 
mation, we obtain the generalized Gross-Pitaevskii (GP) equation, 

h 2 d 2 



uj 2 z 2 + K P + e« 



2m dz 2 



+ E 9aa'/3/3' 



(6) 



The Popov approximation neglects the anomalous correlation ypipj jl_5[] . From the numerical 
solutions of the GP equation ([6]) using the trap frequencies relevant to the experiment Q], 
we find that |$ Q | 2 /|$o| 2 ^ 10~ 6 (a ^ 0), where we have denoted the lowest radial state 
as a = = (0,0), and thus the contribution from higher radial modes to the condensate 
wavefunction is negligible small. For this reason, we will henceforth approximate $ Q ~ $# Q ,o- 
Taking the usual Bogoliubov transformations for the noncondensate, 



(?) 



we obtain the coupled Bogoliubov equations, 



L a Uja + E a ' ( 2 9a' n + 9aa' (3p>npp>) U ja > ~ g^UQVj^ 



(8) 



L a Vja + Ea' (^9a' n + 9 W W^w) V 3<*' ~ 9a> n Q u ja' 



(9) 



where we have introduce the operator 




(10) 



We have also introduce the simplified notations n (z) = |<fr(z) | 2 , g%, = g aa 'ooi nppi = 



As noted above, we only include the lowest mode (a = 0) in the condensation wavefunction 
$. Sums over the repeated indices /3,/3' are implied in Eqs. ([8]), and ([9]). These equations 
define the quasi-particle excitation energies Ej and the quasi-particle amplitudes Uj a and 
Uj a . The orthonormality of the quasi-particle amplitudes is specified by the relation 



Using the solutions of Eqs. ([8]) and Q, one can obtain the noncondensate density from 
n — E Q n aa , where 



with N(Ej) = 1/ [expipEj) - 1}. 

Solving the coupled equations (jB]), (|S]) and flU]), we self-consistently determine the exci- 
tations spectrum Ej and the condensate fraction at finite temperatures. Our calculation 
procedure is summarized as follows. Eq. ([6]) is first solved self-consistently for \x and $ ne- 
glecting the interaction terms. Once $ is known, Ej, Uj a and Ujp are obtained from Eqs. dHJ 
and (jUJ) with hppt set to zero. This is inserted into Eq. (jSJ) and the process is repeated until 
convergence is reached. At each step, we define the normalization of the condensate wave- 
functions by/ dz |$(z)| 2 = N—N, where N = J dz n(z) is the total number of noncondensate 
atoms. 

Throughout this paper we use the following parameters of the experiment of Ref. {4]: 
m (87 R b)=i.44xl0- 25 kg, u z /2tt = 9.0 Hz, u ± /2tc = 92 Hz, scattering length a = 5.82 nm 




(11) 



*M = Z ( U 3^j/3 + V 3<± V 3P) N ( E j) + V j*Vj(3 



(12) 
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and the wavelength of the optical lattice A=795 nm. We fix the total number of atoms as 
N = 4 x 10 5 . In Fig. [H we plot the condensate fraction N c / N as a function of the temperature 
for various values of the lattice depth s by solving GP equation Eq. ^ and Bogolibov 
equations Eqs. (jSj) and 0. It can be seen that each line falls to zero at approximately 
T =140 nK, which is close to the semiclassical prediction of the BEC transition temperature 
of an ideal Bose gas in a 3D harmonic trap T c °= QMhiojlu^N 1 / 3 /k B = 141 nK In 
contrast, T c of an ideal Bose gas in a ID harmonic trap is T c ~ %uj z ^^^ =13.4 /iK jl?! ]. 
This crearly shows that one must explicitly take into account the excitations in the radial 
direction in order to obtain the correct thermodynamic behavior at finite temperatures 
[14J]. The transition temperature deceases with increasing lattice depth, but there is no 
significant change in the temperature dependence of the condensate fraction, as long as the 
lattice potential is not so deep, i. e., s < 2, 
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FIG. 1: The condensate fraction N c /N as a function of the temperature. "s=0" represents the 
ideal Bose gas result without a lattice potential N c /N = 1 — (T/T,?) 3 , where T® =141 nK. 



We turn to the detailed structure of the excitation spectrum Ej. From Eqs. © and 
([9]), one sees that the different radial modes are coupled due to interactions. However the 
coupling is so small that there are still distinct branches corresponding the radial modes. 
We thus label each branch with the index a. For example, the lowest branch (a = 0) can be 
identified with the branch corresponding to the lowest radial mode. In Fig [21 we plot the 
excitation spectrum. Frequency of the condensate collective mode is related the excitation 
energy through Ej = hcuj. In particular, the dipole mode frequency can be identified with 
the lowest frequency. In Fig. [3l we plot the dipole-mode frequency as a function of the 
lattice depth with a fixed temperature T= 40nK. One can see the decrease of the dipole- 
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FIG. 2: The Bogoiubov excitation energy spectrum for s=1.2, where i is the energy label for each 
branch 



mode frequency with increasing lattice depth. The same behavior was also observed in 
the experiment jsj]. The negative energy shift also found for other excitation modes. The 
decrease of T c can be attributed to these negative shifts. 
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FIG. 3: The dipole mode frequency as a function of the lattice depth s with a fixed temperature 
T = 40nK. 



III. LANDAU DAMPING RATE OF DIPOLE MODE 

In this Section, we calculate Landau damping of the dipole mode. Landau damping is 
the dominant damping mechanism for lowdying collective modes in trapped Bose-condensed 
gases in the collisionless regime at finite temperatures ja|. Landau damping originates from 
the coupling between single-particle excitations and the collective oscillations. Damping 
occurs because the thermal bath of the elementary excitations can absorb quanta of the 
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lective oscillations. A general expression for the Landau damping rate is given in Refs 
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.a 



17| . In our quasi-lD model, Landau damping rate can be expressed in terms of the 



quasi-particle excitation energies E~ and quasi-particle amplitudes Uj a and Vj a calculated in 



the previous Section: 



7 l = 4vr £ Wo0 £ \A«« | 2 [iV(£ 4 ) - N^E^Qiw + E t - E 3 ) 

a,a' i^j 



(13) 



where uo is the eigenfrequency of a condensate dipole oscillation. The matrix element A 
is defined by, 



act 

ij 



-ViO [u* a Uj a > + v* a v ja > - u* a v jal }} 



(14) 



where u w and t>i are the quasi-particle amplitudes corresponding the dipole mode. One 
difficulty in calculating the damping rate in Eq. (1131) using the discrete excitation spectrum 
Ej is that it involves the energy-conserving delta functions. This difficulty can be overcome 
by replacing each delta function by a function with a finite width. In this paper, we use the 
following replacement: 



5 (hu + Ei- Ej) -► — 6 ( A-hu + Ei 
2 A v 



E, 



(15) 



The width factor A is somewhat arbitrary, and the result for 7l will vary with A. In Fig. HI 
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FIG. 4: The variation of the Landau damping rate 7l (arbitrary unit) with the width factor A 
for lattice depth s = 1.0 and temperature T = 120 nK. To capture the A dependence of 7l while 
saving computational time, we take the matrix element Af- 1 ' as a constant in this figure. 



we plot a typical behavior of the A dependence of 7l- One can see that the variation of 
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7l is very weak when A/hu z lies between lxlO -4 and 3 xlO~ 4 . The same behaviors are 
also found for other temperatures and lattice heights. We thus take A/hu z = 2.5 x 10~ 4 in 
calculating the damping rate 7l- 
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FIG. 5: The damping rate as a function of lattice depth s with a fixed temperature T = 120nK. 
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FIG. 6: The damping rate as a function of temperature with a fixed lattice depth s = 1.0 

In Fig. [5], we plot the damping rate as a function of lattice depth s with a fixed temperature 
T = 120nK calculated from Eq. ([TBI . We see that the magnitude of the damping rate is 
~ Is -1 , which is in reasonable agreement with the experimental data [J]. We also see 
that the damping rate increases with increasing lattice depth, which is consistent with the 
experimental result |4J]. This increase of the damping rate can be attributed to the increase 
of the number of the elementary excitations because of a reduction of the excitation energy 
with increasing lattice depth. 

We next investigate the temperature dependence of the damping rate. In Fig. [61 we plot 
the damping rate as a function of temperature with a fixed lattice depth s = 1.0. We find that 
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damping rate decreases rapidly with decreasing temperature. This is due to the reduction 
of thermal excitations with decreasing temperature. Because of this strong temperature 
dependence of the damping rate, it is very difficult to make a quantitative comparison with 
the experimental data without precise knowledge of the experimental temperature. 

Here we comment on the effect of adiabatic loading of an optical lattice on the temperature 
of a gas. The usual path for preparing for condensed Bose gases in an optical lattice consists 
of first forming a ultracold bosons in a weak magnetic trap, to which a ID lattice potential 
is adiabatically applied by slowly rating up the light field intensity. In this case, the initial 



and fina 
loading 



temperature are not usually equal since the energy spectrum changes during lattice 



161 ] . We used the entropy-temperature curves to consider the effect of the adiabatic 



loading into a lattice. For the initial (s=0) temperature T=120nK, the temperature shift in 
the final state s = 1.6 is only about 2nK. Thus, one can ignore the effect of the adiabatic 
loading in shallow lattice regime. 

We note that the radial excitations are important in Landau damping. If we calculated 
the the damping rate ignoring the radial thermal excitations, the damping rate would be 
order of magnitude smaller. This means that the radial thermal excitations make significant 
contributions to the Landau damping. 



IV. CONCLUSION 



In this paper, we studied the Landau damping of dipole oscillations of Bose-condensed 
gases in a combined potential of highly-elongated harmonic trap and ID optical lattice, with 
explicitly including the effect of the radial excitations. While we treated the condensate 
wavefunction only with the lowest radial mode, we took into account the radial excitations 
for thermal cloud. 

First, we studied equilibrium properties of Bose-condensed gases in a combined harmonic 
trap and ID optical lattice potentials. We have presented a detailed calculation of the 
condensate fraction in a ID optical lattice at finite temperatures. We find the negative shift 
of the excitation energies with increasing lattice depth. We obtain the dipole-mode frequency 
as a function of the lattice depth. The negative shift of the dipole-mode frequency was also 
observed in the experiment. 

Second, we calculated Landau damping rate of dipole-modes with varying lattice depth 



10 



and temperature. Result for the damping rate is consistent with experimented data. There- 
fore, the experimentally observed damping can be understood as Landau damping. We also 
showed that the radial thermal excitations are important in both equilibrium condensate 
fractions and Landau damping rate. 
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